gusucode.com > MATLAB与控制系统仿真实践 源代码 > MATLAB与控制系统仿真实践 源代码/第十四章/eg14_6.m

    num=100;
den=[1 8 0];
w=logspace(-1,2,100);                     %指定频率范围
nichols(num,den,w);                        %绘制系统的Nichols曲线
ngrid;

%%%%%%%%%%%%%%%%%%%%%%%%
num=100;
den=[1 8 0];
w=logspace(-1,2,100);
[mag,phase]=nichols(num,den,w);                      %返回Nichols曲线参数
magdB=20*log10(mag);                                      %幅值单位转换
subplot(2,1,1)
semilogx(w,magdB);                                           %使用Nichols曲线参数绘制幅频特性
title('系统幅频特性曲线');
subplot(2,1,2)
semilogx(w,phase);                                              %使用Nichols曲线参数绘制相频特性
title('系统相频特性曲线');
figure(2);
bode(num,den,w)                                                  %求取系统Bode图
title('系统Bode图');